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Abstract 

A method is presented for approximating the effective conductivity 
of composite media with thin interphase regions. The methods pre- 
sented are computationally efficient in the sense that the fields need 
to be computed only in a reference composite in which the interphases 
have been replaced by perfect interfaces. The results apply whether 
any two phases of the composite are separated by a single interphase 
or multiple interphases, whether the conductivities of the composite 
phases are isotropic or anisotropic, and whether the thickness of an 
interphase is uniform or varies as a function of position. It is assumed 
that the conductivities of the interphase materials have intermediate 
values as opposed to very high or very low conductivities. 
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1 Introduction 



Composite media with thin interphases between adjacent phases are used 
widely in engineering, biomedical, and technological applications. Examples 
include coated fiber or coated particle reinforced materials and composites 



with adhesively bonded joints, see for instance Kim and Mai (1998); Banea 



and Da Silva (2009); Ramakrishna et al. (2001). In this work we consider 



the problem of computing the effective conductivity in a composite with thin 
interphase regions between some of the composite's constituents, see Figure[T] 
The fact that the composite has a thin interphase (or interphases) makes the 




Figure 1: (a) A composite in which media 2 and 3 form an interphase region 
between media 1 and 4. (b) A reference composite in which the interphases 
of the composite in (a) have been replaced by an interface. 



numerical computations of the fields and the effective properties a difficult 
task. Therefore, approximate models of thin interphase composites that do 
not require solving for the fields inside the interphase regions are desirable, 



Benveniste (2012). 



There has been an extensive study of interface models, in which an inter- 
phase region between two media is replaced by an interface, allowing direct 
contact of the media, along with appropriate interface conditions. This inter- 
face is usually referred to as an imperfect interface when either the potential 
or the normal flux component has a jump discontinuity across the inter- 
face. The effective conductivity in composites with highly conducting or 



poorly conducting imperfect interfaces has been studied by Benveniste and 



Miloh (1986); Benveniste (1987); Torquato and Rintoul (1995); 



Vernescu (1996); Lipton (1997); Miloh and Benveniste (1999); 



Lipton and 



Lipton and 



Talbot (2001); Le-Quang et al. (2010), among others. Numerical methods 



for computing the effective conductivity in composites with highly conduct- 
ing or poorly conducting imperfect interfaces have been recently developed 
by |Yvonnet et al.| ( [20081 

Although the majority of studies in the literature of thin interphases fo- 
cus on the case when there is a high contrast between the conductivities of 
the interphase and its adjacent media, there have been some works that deal 
with intermediate values of conductivities without the high-contrast assump- 
tion. Notable examples of such studies include the work of Hashin (2001 ) and 



the works of Benveniste (2006a 



(2012). In the work of Hashin 



b|£ Benveniste and Baum (2007): 



Benveniste 



(2001) an interface model is introduced, in 
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which the interphase's conductivity is arbitrary ranging from zero to infinity, 
and the effect of the interphase on the effective conductivity is discussed. 
However, this work is restricted to having isotropic and homogeneous phases 
and the thickness of the interphase is assumed to be constant. In the above 
studies by Benveniste et. al. two models of thin interphase materials are 
developed of which one is an interface model with the interphase replaced by 
an interface. In the second model the geometry of the thin interphase is left 
intact and is characterized by conditions in which the fields are evaluated in 
the adjacent media at both sides of the interphase and do not involve the 
fields within the interphase. Some of the interesting features of these models 
are that the interphase can have anisotropic and inhomogeneous conductiv- 
ities. However, it is assumed that the thickness of the interphase region is 
constant. It is also noted that these studies do not discuss the effect of the 
thin interphase on the effective conductivity. For a more comprehensive re- 
view of the topic (on the modeling of thin interphases), we refer the reader 



to the works given in 



Benveniste 


(2006a 


); 


Gu and He (2011) 



In this work we present a novel method that provides an approximation 
to the effective conductivity and combines different desirable features. These 
include that the method is computationally inexpensive in the sense that 
it does not involve computing the fields within the interphase region. Ad- 
ditionally, the method applies for homogeneous as well as inhomogeneous 
interphases, that is whether a single interphase separates two materials or 
multiple (possibly infinitely many) interphases. Moreover, the conductivities 
of the composite components, including the interphases, can be isotropic or 
anisotropic. Furthermore, the geometry of the interphase is arbitrary and its 
thickness h is not assumed to be constant but can vary slowly, on a length 
scale large compared to h. However, the total thickness of the interphase re- 
gion must be sufficiently small in order for our approximation of the effective 
conductivity to hold. It is also assumed that the interphase conductivities 
have intermediate values as opposed to very high or very low conductivities. 
Numerical results for a particular example of a thin-interphase composite 
show that our approximation agrees very well with the exact effective con- 
ductivity when the interphase is thin and has intermediate values of conduc- 
tivity. 

Our method is based on estimating the change in the effective conduc- 
tivity as the total thickness of the interphases vanishes. Thus we compare 
the effective conductivity of the composite to a reference composite in which 
the interphase materials in the original composite have been removed and 
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replaced by a perfect interface, with the usual conditions of continuity of po- 
tential and continuity of flux across this interface . Computing the effective 
conductivity and the fields inside the reference composite is an easier problem 
that may, for example, be solved using integral equations. We refer to this as 
the unperturbed problem. The method is based on the following estimate, 
which holds when the total thickness of the interphases is sufficiently small, 

<r* « a* = cr* + 5(7*. (1.1) 

Here cr* is the effective conductivity of the composite with thin interphases, 
cr* is our approximation, cr*, is the effective conductivity of the identical 
composite but with the interphase materials removed and replaced by one 
material of the two neighboring phases, and 5cr* is the change in the effective 
conductivity due to inserting the interphase materials. The formulas that we 
present and derive in Section [3] show that in order to compute 5cr* we need to 
compute the fields only at material interfaces inside the unperturbed problem 
(or reference composite). 

This article is organized as follows. In Section |2j a review is given for com- 
puting the change in the effective conductivity due to translating an interface. 
In Section |3j a formula for computing the change in the effective conductivity 
due to inserting an interphase region is derived. Finally, Section [4] provides 
numerical results for the effective conductivity of the doubly coated sphere 
assemblage and shows a comparison between our approximation, the exact 
result, and two other known approximations. 



2 Review 



In this section we review how the effective conductivity of a composite 
changes due to a shift of a phase boundary. The review here is based on 



the treatment given in Section 16.6 of Milton (2002). For completeness of 



the presentation, we will include a derivation of Milton's interface-shift for- 



mula (2.12), in the case of the conductivity of composites. 



Consider a periodic composite with a period cell given by Q C 



Let 



r = be an interface between two phases, where r\ is a continuous parameter 
such that the position of the interface r v changes when r\ changes. We denote 
by cr,, the conductivity inside the composite and assume that cr^x) = cr* 
for x inside the i-th phase, where cr 1 is a symmetric matrix representing an 



4 



- side f n 




Figure 2: Perturbation of an interface separating two phases. The distance 
between the interfaces T v and T v+e is approximated, to the first order in e, 
by e l7r?(x)|. 

anisotropic conductivity. The electric and current fields are denoted by 
and J v , respectively, and satisfy the differential constraints, 

V-J r , = 0, E v = Vu v , J v = cr v E v , (E v ) = E , (2.1) 

where u v is the electric potential and E is the applied electric field which 
is chosen to be independent of 77. Here, and in the rest of this article, the 
brackets (•) denote taking the average over fl, that is, 




We assume that the potential u v and the flux n ■ 3 V are continuous across 
the interface T v , where n is the normal to the interface. These continuity 
conditions can be explicitly written as 

u+(x) = u~(x), n(x) ■ J+(x) = n(x) ■ J~(x), for all x G T v , (2.2) 

where the plus and minus superscripts denote the two sides of the interface 
with the convention that the normal is directed outward from the + side 
pointing toward the — side, see Figure [2} 

For a small perturbation in rj, rj — > rj + e, the translated interface r^ +e 
can be approximated to the first order in e by 

T 'ri,e = { X V = x + e7 r? (x)n(x), x G T v }, (2.3) 
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where 7^ is a scalar such that the quantity |e7^(x)| represents the distance 
traveled by the point x e T v along the direction of n(x). We note here that 
the fact that the potential u v+e is continuous across r,j +£ implies that, to the 
first order in e, 



u 



x + e7 r) (x)n(x)) 



u 



x + e7 ?? (x)n(x)). 



(2.4) 



Let <t* be the effective conductivity of the composite. Then the derivative 
of the energy with respect to rj is given by 



d 
drj 



(Eq • c*E ) 



d 
dr] 
2 

W\ 



dE v 

n\ r„ dr] 



+ \k I -*K-* +E J 

I I J T v 



El cr E 



(2.5) 



where the first integral represents the change in the energy due to variations 



in the field 



dE, 

dr) 



treating the interface as fixed, and the second integral 



represents the change in the energy due to translating the interface, treating 
the field E^ as fixed outside the neighborhood of the interface. 
Integrating by parts, the first integral in (2.5) becomes 

duZ 



dE v 

dr] 



9ut , , N 
V n • J+ - 



dr] 

dul 



v 
duZ 



dr] dr] 



fn • J 



(2.6) 



where in the last equality we used the continuity of the flux at the interface 
which is given by (2.2). By expanding (2.4) and equating the terms of order 
e, we obtain 



duZ) 

— + 7 ,n- Vu v 



du 



dr] 



— + 7^11 • Vu . 



(2.7) 



By combining (2.7) and (2.6), we see that 

/ 7^n- (Vm+- Vm^) (n- J; 



r dE 

Jn\ r v dr] 



/ In (V< - Vu") • (n ® n) J; 



(2.8) 



Let Wi(x), w 2 (x), . . . , W7 rf _x(x), and n(x) be an orthonormal basis at each 
point x e T v . Then since the potential is continuous across the interface, 
it follows that the tangential derivative of the potential is also continuous 
across the interface and hence 



W; 



(Vm+-V«-)=0, fori = l,2,...,d-l. 



Using this fact, (2.8) can be written as 

r dE 

Jn\ r n 



j 7, ( v < - v v ) • (E 

/ In (V< - Vu~) ■ J + 
j 7 „(E--E+) 



Wi ® Wi + n g) n) J; 



cr+Ej 



(2.9) 



Using the same arguments above, only replacing n • J| by n ■ J in (2.6), 
gives 



n\ r, 



cr 



, = / 7, (E- - E+) 
•/r n 



c-E- 



(2.10) 



By substituting (2.9) and (2.10) back into (2.5), we obtain 



da* 
E ■ — — E 

drj 



d_ 
drj 
1 



(E • a*E ) 
^ jf l v i^-cT + E+-E+- ( T-E; ri 



E v ■ 3 V \. 



(2.11) 



The following theorem summarizes the results of this section. 

Theorem 2.1 QMilton| 02002) ). A ssume that an interface T is translated by 
a distance e 7 (-) in the direction normal to the interface, keeping the applied 
electric field (E) fixed. Then the resultant change in the energy, to the first 
order in e, is given by 



(E) • 5a* (E) » p J e 7 (x) [E~(x) ■ J + (x) - E+(x) ■ J~(x)] dx, 

where 5a* denotes the change in the effective conductivity a* . 



(2.12) 
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3 Change in the effective conductivity of thin- 
interphase composites 

Consider a multi-phase composite consisting of a finite number of homoge- 
neous constituents. Let cr° and cr 1 be the anisotropic conductivities of two 
of these constituents such that between them lies a set of K — 1 thin inter- 
phases, whose total thickness is h, having anisotropic conductivities given 
by cr 2 , cr 3 , . . . , a K . We analyze the change in the effective conductivity of 
this composite as the interphases total thickness vanishes (h — » 0). In the 
limit, the interphases are replaced by an interface between the 0-phase and 
the 1-phase of this composite. 



In Sections |3.2| and |3.3| we present and derive formulas for this change in 
the effective conductivity for the case of a single thin interphase and for the 
case of two thin interphases, respectively. The general case of a finite set of 



thin interphases is discussed in Subsection 3.4, Section 3.5 focuses on the 



case when the interphase conductivity is continuously varying. 
3.1 Preliminaries 

From Theorem 2.1 translating an interface T by a distance h in the direction 
normal to the interface, keeping (E) fixed, results in a change in the energy, 
to the first order in h, given by 

(E) ■ 5a* (E) « / h [E~ • J + — E + • J~] . (3.1) 
I" I Jr 



Note that the integration variable in (3.1) has been suppressed for concise- 



ness. This result can be generalized to include the case of simultaneously 
shifting more than one interface. Specifically, if a finite number of interfaces 
Ti, . . . , Ttv are shifted by distances t\, . . . , ijv, in the directions ni, . . . , n^v, re- 
spectively, where rij is the direction normal to interface Tj, then the resultant 
change in the energy is given by 

- 1 f 

(E) • 5a* (E)*J2m [ET • J t ~ E t ■ 37] ■ (3-2) 
i=i JTi 

For a field F, we denote by F t and F n the tangential and normal compo- 
nents of the field, respectively. Using this notation, we present the following 



lemma which provides an equivalent yet useful way for writing (3.1). 
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Lemma 3.1. Let cr + and cr be the conductivities of phase + and phase — 
as described above. Then the change in the effective conductivity is given by 



(E)-5<r*(E) « j^r f h [E~ • J + - E + • J~] 

= j^|^[(o- + - --)E t + -E+-(( ( 



,-1 



(3.3) 
(3.4) 



Proof. We expand the integrand by writing the fields in terms of their tan- 
gential and normal components to obtain 



E J 



(E,-+E-)-(J+ + J+) 

E r • Cr+E ^ + E iT • J n + E n ' J t" + E n ' J n 



(3.5) 



Using the facts that E t = Ef, J n = J+, and E n = (cr ) 1 J n , (3.5) becomes 
E--J + = cr+E+ • E+ + (cr-)~ l J+ • J+ (3.6) 

Similar arguments give 

E+-J- = cr-E+ ■ E+ + (cr+)- 1 J+ • J+ (3.7) 



Taking the difference between (3.6) and (3.7) and substituting the result in 
(3.3) gives (3.4), completing the proof. □ 



Remark. 



Since E t = E^,J n = J+, (3.4) can equivalently be written as 



(E) • 5cr*(E) 



1 



J h [(cr + - 0-)E- ■ ET - ((CT+)- 1 - (cr-)" 1 ) J n ■ J~] 



3.2 Single interphase 

Consider a composite in which two phases, with conductivities cr° and cr 1 , are 
separated by an interphase of conductivity cr 2 and thickness h, see Figure [3] 
We note that the thickness h is not assumed to be constant but can vary 
slowly along the interphase region. Let F± and T2 denote the interfaces 
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Figure 3: A single-phase interphase region. In the limit as I\ — > T 2 , the 
interphase vanishes and is replaced by the interface T 2 which now separates 
phase-0 and phase- 1. The figure is schematic in that the thickness h is not 
assumed to be constant but can vary on a length scale large compared to h. 



between phase-1 and phase-2, and phase-2 and phase-0, respectively. We 
may assume, without loss of generality, that the directions of the normal 
vectors ri! and n 2 are given according to Figure [3] We also assume that 
the interphase geometry is locally laminar, which in particular implies that 
ni approaches n 2 , as h approaches zero. Then by applying Lemma 3^ we 
find that the change in the effective conductivity due to translating Ti by a 
distance h in the direction of rii is 

(E)-5o-*(E) w / h [E 2 ■ J 1 — E 1 ■ J 2 ] 

(3.8) 



3.3 Two interphases 

Consider a composite in which two phases, with conductivities cr° and cr 1 , 
are separated by two interphases of conductivities cr 2 and cr 3 , and thickness 
h 2 h and h 3 h, where h 2 + h s = 1, see Figure [4} We note that the thicknesses 
h, h 2 h and h 3 h are not assumed to be constant but can vary slowly along the 
interphase region. Let r 1; T 2 , and T 3 denote the interfaces between phase-1 
and phase-2, phase-2 and phase-3, and phase-3 and phase-0, respectively. 
We may assume, without loss of generality, that the directions of the normal 
vectors ni, n 2 , and n 3 are given according to Figure |4| We assume further 
that the interfaces Ti and T 2 are translated simultaneously by distances 
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h = hqh + h 3 h and h^h in the directions xii and n 2 , respectively. We again 
assume that the interphase geometry is locally laminar, which implies that 
ni and ri2 approach 113. Then we have the following result. 

Theorem 3.1. The change in the effective conductivity 5cr* due to the change 
in the composite as described above is given through the formula 

<E)-5<7*(E> « ^-J^h[{cr l -{h 2 a- 2 + h^))E)-E\ 

- ii* 1 )- 1 ~ {hii* 2 )- 1 + hi* 3 )- 1 )) J* ■ Ji]. (3.9) 



Proof. By using Lemma 3.1 and (3.2) we obtain that the resultant change in 
the effective conductivity is given by 

(E)-8a*(E) « j^r f h [E 2 ■ J 1 — E 1 • J 2 ] + / hh 3 [E 3 • J 2 - E 2 • J 3 ] 
"I Jti |"| </r 2 



1 



+ M /r " ^ ' E? " {{(T2yl ~ 31 ' 3l 

(3.10) 

By taking the limit as h — y 0, we have that T 2 ~ r 1; E 2 w Ej, and J 2 « J„. 
Therefore, we can write (3.10) as 

(E) ■ 5a* (E) « J_ ^ /, [ ((T i _ CT 2 + ^2 _ ^ 3(T 3 )E i . E i 



2\-l 



(3.11) 



From (3.11) and by using the fact that 1 — /13 = /12, equation (3.9) follows. □ 



3.4 Multiple interphases 

Consider a composite in which two phases, with conductivities cr° and cr 1 , 
are separated by K interphases of conductivities cr 2 , cr 3 , . . . , cr K , and thick- 

We note that the 
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ness h 2 h, h 3 h, . . . , where /;.£; = 1, see Figure 



k=2 



h A h 



i 




r, 



n 



n 



Figure 4: A two-phase interphase region. In the limit as I\ — > T 3 and 
^2 ^ T 3 , the interphase region vanishes and is replaced by the interface T 3 
which now separates phase-0 and phase-1. The figure is schematic in that 
the thicknesses h, hih and h%h are not assumed to be constant but can vary 
on a length scale large compared to these thicknesses. 

thicknesses h, hih, h^h, . . . , h^h are not assumed to be constant but can vary 
slowly along the interphase region. Let Ti, r 2 , . . . , Tk denote the interfaces 
between phase-1 and phase-2, phase-2 and phase-3, . . . , and phase- K and 
phase-0, respectively. We may assume, without loss of generality, that the 
directions of the normal vectors iij, n 2 , . . . , are given according to Fig- 
ure |5} We assume further that the interfaces T ly . . . ,T K _i are translated 
simultaneously such that each F k , k — 1, . . . , K — 1, is translated by a dis- 

K 

tance h hi in the direction normal to the interface T^. Then we have 

i=k+l 

the following result. 

Theorem 3.2. The change in the effective conductivity 5cr* is given through 
the formula 




(3.12) 
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Figure 5: A multi-phase interphase region. In the limit as Ti — > Fx, 
^2 — > F K , . . . , r^_i — > T K , the interphase region vanishes and is replaced 
by the interface IV which now separates phase-0 and phase- 1. The figure is 
schematic in that the thicknesses h, Ii2h, h^h, . . . , hxh are not assumed to be 
constant but can vary on a length scale large compared to these thicknesses. 



Proof. From (3.2) we see that the resultant change in the energy is 

K-l (K-l \ 

(E)-<5<r*(E) « J2jo\ h hi )[ Ek+1 ■ Jh ~ Efc • Jk+1 ] ■ 

k=l I I J T h \i=k+l J 



(3.13) 



Then, by using Lemma 3.1, (3.13) can be rewritten as 

K-l 



k=l 1 1 Jl k \i=k+l 



hi I (a k - a k+1 )E k t ■ E k t 



"((« 



f^k+l\-l\ -ik -rk 



(3.14) 



In the limit as h — >■ 0, we have that T k ~ fi, « Ej, and J„ « J„, for 
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k = 2, 3, . . . , K — 1. Therefore, we can rewrite (3.14) as 
(E)-5<r*(E) ? 



4 / ME ( E (- fc - - fc+1 ) E * • E * 

I I fc=i \i=fc+i / 



K-1 / 



-^E E * ((^)- 1 -(^ +1 r 1 )^-^: 



k=i \i=k+i 



(3.15) 



A" 



By expanding the sums and using the fact that hk — 1, we find that 

fc=2 

E E M<>- fc -0- fc+1 ) = ^-E^ (3-16) 

k=l i=k+l k=2 

EEMW -1 -^ 1 ) -1 ) = (O^-EWr 1 . (3.17) 



fc=l i=fc+l 



fc=2 



The result follows from substituting (3.16) and (3.17) in (3.15). 



□ 



3.5 Graded interface 



By taking the limit of (3.12) as K — > oo, one obtains the following result for 



an interface with graded conductivity 



(E) • 8a* (E) 



1 



ri L 



h a 1 



a(z) dz ) Ej • Ej 



(cr)- 1 ^) dz ) Ji -4 



. (3.18) 



It is worth noting that this formula applies even if a has a finite jump at 
z = or z = h. 



4 Numerical Results 

In this section, we consider a specific example of a composite with an in- 
terphase; namely the doubly coated sphere assemblage for a three-phase 
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Figure 6: A doubly coated sphere from the doubly coated sphere assemblage. 



material, see for instance Milton (2002). The reference composite is the 



well-known singly coated sphere assemblage of Hashin and Shtrikman Hashin 



(1962); Hashin and Shtrikman (1962). The goal is to provide a numerical 



comparison between our method for approximating the effective conductivity 
and the exact effective conductivity of the doubly coated sphere assemblage. 
The numerical comparison also includes two other well-known approxima- 
tions for the effective conductivity for the cases of highly conducting inter- 
phase and poorly conducting interphase. 

For convenience, we assume that the phases of the doubly coated sphere 
assemblage have isotropic conductivities and hence, for this composite, the 
effective conductivity is also isotropic. Therefore, in this section, we represent 
conductivities by scalar quantities and denote by a*, 5cr*, cr°, and c^, for 
i = 1,2, 3, the effective conductivity, the change in the effective conductivity, 
the effective conductivity of the reference composite, and the conductivity of 
the i-th phase, respectively. 



4.1 The doubly coated sphere assemblage 

Consider an arbitrary doubly coated sphere in the doubly coated sphere as- 
semblage in which the three concentric spheres have radii given by n , r?, , 
with 7*1 < r 2 < r 3 , see Figure [6] Assume that the core is occupied by a ma- 
terial with conductivity o"i, the interphase is occupied by a second material 



15 



cr, = 1, a = 10, 




Figure 7: Log plot of the effective conductivity for the doubly coated sphere 
assemblage as a function of the interphase conductivity a%. The interphase 
volume fraction is 82 = 0.1. The core volume fraction is Q\ = 3 3 /4 3 w 0.422, 
which is the same in both the original and reference configurations. The 
outer coating and the core have conductivities <j\ = 1, cx 3 = 10, respectively. 
The interphase conductivity values are in the range 02 G [10~ 2 , 10 2 ] with a 
zoomed- in view in the middle for conductivity values a 2 G [1, 10]. 



with isotropic conductivity cr 2 , and the outer-most coating is occupied by a 
third material with isotropic conductivity 03. With this notation, the effec- 
tive conductivity of the doubly coated sphere assemblage can be computed 
explicitly 



see for example, 
( 2002[ )) and is given by 



Schulgasser (1977) and Section 7.2 of Milton 



3<r 3 (l-03 



(4.1) 



cr 3 - cr 2 



9 2 - 



3(7901 

3<r 2 (l-fl 3 

02 - 01 



16 



where 




are the volume fractions of phase- 1, phase-2, and phase-3, respectively. 

It is important to note that the thickness of the interphase, phase-2, 
varies from one doubly coated sphere to the next doubly coated sphere, in 
proportion to the relative sizes of the doubly coated spheres, so that the 
doubly coated spheres are rescaled copies of each other. 



4.2 The approximate effective conductivity 



In the doubly coated sphere assemblage described in the previous section, 
let ri be fixed and r 2 = (1 + t)r\ for some t > 0. Then the thickness of 
the interphase is h = r 2 — r\ = tr\. By letting t — > 0, the thickness of the 
interphase h 0, and hence we can compute the change in the effective 

However, in this case it is more convenient 
5a* da* 

dh 



conductivity 5a* using (3.8). 



and that a* is given 



h=0 



to compute 8a* from noticing that lim 

/i^o h 

explicitly by (4.1). Thus, by using the fact that h = r 2 — r\, substituting 
(4.2) in (4.1), and then differentiating with respect to h, one obtains 

da* 



8a*(h) 



h 



dh h=0 
9a 3 9i(a 3 



^)((o r i) 2 + 2a 2 a 3 ) 



h 1 (4.3) 

na 2 ((a 3 - ai)9 1 + a 1 + 2a 3 ) 2 

On the other hand, the effective conductivity of the reference configuration 
(the unperturbed problem) cr° can be calculated through evaluating a* at 
h = 0, which after simplification gives 

3a 3 ^i 



cr 3 + 



1 



3cr 3 
a 3 - oi 



(4.4) 
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Figure 8: Log plot of the effective conductivity for the doubly coated sphere 
assemblage as a function of the interphase conductivity 02. The interphase 
volume fraction is 62 = 0.01. The core volume fraction is 9\ = 3 3 /4 3 w 0.422, 
which is the same in both the original and reference configurations. The 
outer coating and the core have conductivities <j\ — 1, <t 3 = 10, respectively. 
The interphase conductivity values are in the range 02 G [10~ 2 , 10 3 ] with a 
zoomed- in view in the middle for conductivity values 02 G [1, 10]. 



By substituting (4.3) and (4.4) in the scalar version of (1.1 ), one obtains the 



following formula for our approximation of the effective conductivity for this 
composite 







.(h) 



cr*(h) 
0-3 + - 



3a 3 #i 



01 



3^3 

0-3 - 01 



+ h 



-9cx36>i(cx3 - a 2 )((ai) 2 + 20203) 



ria 2 ((o-3 - cri)0i + en + 20 3 ) 5 



(4.5) 
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a 1 = 10, a 3 = i, e = 0.1 




10° 10' 

U 2 



Figure 9: Log plot of the effective conductivity for the doubly coated sphere 
assemblage as a function of the interphase conductivity a-i. The interphase 
volume fraction is 82 = 0.1. The core volume fraction is Q\ = 3 3 /4 3 w 0.422, 
which is the same in both the original and reference configurations. The 
outer coating and the core have conductivities a\ = 10, cr 3 = 1, respectively. 



4.3 Approximations of the effective conductivity for 
high-contrast cases 



We consider the doubly coated sphere assemblage described in Section [41] for 
the cases in which the interphase is thin and is either highly conducting or 
poorly conducting. The high-conduction case corresponds to assuming that 
#2 — and cr 2 — > 00 in such a way that the product #202 remains constant. 



By taking these limits in (4.1), we obtain that the effective conductivity in 
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a l = 10, a } = i, 6 = 0.01 



2.65 



2.55 



b 2.5 




2.45 



2.35 



Figure 10: Log plot of the effective conductivity for the doubly coated sphere 
assemblage as a function of the interphase conductivity o^- The interphase 
volume fraction is 62 = 0.01. The core volume fraction is 9\ = 3 3 /4 3 w 0.422, 
which is the same in both the original and reference configurations. The 
outer coating and the core have conductivities a\ = 10, cr 3 = 1, respectively 



this case is given by 

_ „_ a. 3 (!- 

3o- 3 



■f = *.+ • (4-6) 



03- 



2^ 2( T 2 

(T3 - (71 



3(1 - 0, 



While the poor-conduction case corresponds to assuming that 62 — > and 
o"2 — > in such a way that the ratio U2/O2 remains constant. By taking these 
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a =l,v } = 100, e 2 = 0.1 




Figure 11: Log plot of the effective conductivity for the doubly coated sphere 
assemblage as a function of the interphase conductivity a%. The interphase 
volume fraction is 82 = 0.1. The core volume fraction is Q\ = 3 3 /4 3 w 0.422, 
which is the same in both the original and reference configurations. The 
outer coating and the core have conductivities G\ — 1, 03 = 100, respectively. 



limits in (4.1), we find that the effective conductivity in this case is given by 

3(1 - 3 )a 3 



03 



3a 3 



(4.7) 



0-3 



3 

0"! 



do 



4.4 Numerical comparison of the results 

A numerical comparison given by log plots for the effective conductivity of 
the doubly coated sphere assemblage as a function of the interphase conduc- 
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a, = 100, (T, = I, 



-0.1 



b 3.5 




Figure 12: Log plot of the effective conductivity for the doubly coated sphere 
assemblage as a function of the interphase conductivity 02- The interphase 
volume fraction is 82 = 0.1. The core volume fraction is 8\ = 3 3 /4 3 w 0.422, 
which is the same in both the original and reference configurations. The 
outer coating and the core have conductivities o"i = 100, cr 3 = 1, respectively. 



tivity 02 is provided in Figures [7 12 The five curves displayed in each figure 
correspond to the exact effective conductivity a*, given by (4.1), our approx- 
imation a*, given by (4.5), the highly conducting interphase approximation 



an, given by (4.6), the poorly conducting interphase approximation a^, given 



by (4.7), and the effective conductivity of the reference composite (the zero 
order approximation) cr°, given by (4.4). For a representative doubly coated 
sphere we have taken the outer radius r% = 4 length units and the inner 



radius r\ = 3 length units and thus, by using (4.2), the interphase thickness 
h = r 2 — ri can be computed from the interphase volume fraction 8 2 , which 



we use as a parameter in generating Figures [7 [12 For example, in Figure [7j 
a value of 62 = 0.1 corresponds to h = 0.2204, and a value of 82 = 0.01, as 
in Figure [81 corresponds to h = 0.0235. 
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Figures [7 12 show that there is a good agreement between the approxi- 
mate effective conductivity a* and the exact effective conductivity cr* for the 
doubly coated sphere assemblage. This agreement is most noticeable when 
9 2 , and consequently h, is small, as in Figures [8] and [T0| and for values of 
02 within or near the range [min{ai, CT3}, max{o"i, 03}]. This confirms that 
a* well approximates cr* when the interphase is thin and has intermediate 
conductivity values. 
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